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Employing a nonequilibrium Green's function approach, we examine the effects of long-range 
hole-impurity scattering on spin-Hall current in p-type bulk semiconductors within the framework 
of the self-consistent Born approximation. We find that, contrary to the null effect of short-range 
scattering on spin-Hall current, long-range collisions do produce a nonvanishing contribution to the 
spin-Hall current, which is independent of impurity density in the diffusive regime and relates only to 
hole states near the Fermi surface. The sign of this contribution is opposite to that of the previously 
predicted disorder-independent spin-Hall current, leading to a sign change of the total spin-Hall 
current as hole density varies. Furthermore, we also make clear that the disorder-independent spin- 
Hall effect is a result of an interband polarization directly induced by the dc electric field with 
contributions from all hole states in the Fermi sea. 

b 

PACS numbers: 72.10.-d, 72.25.Dc, 72.25.-b 

I. INTRODUCTION 

Recently, there have been extensive studies of the physics of the spin-orbit (SO) interaction in condensed matter. 
The most intriguing phenomenon induced by SO coupling is the spin-Hall effect (SHE): when a dc electric field is 
applied, the SO interaction may result in a net nonvanishing spin current flow along the transverse direction. The 
SHE is classified into two types according to its origin, an extrinsic spin-orbit Hamiltonian term induced by carrier- 
impurity scattering potentials^ and an intrinsic spin-orbit Hamiltonian term arising from free carrier kineticsi^ The 
intrinsic spin-Hall effect was originally thought to be independent of carrier-impurity scattering. Experimentally, the 
SHE was observed in a n-type bulk semiconductor— and in a two-dimensional (2D) heavy-hole system^ 

However, further studies have indicated that the spin-Hall effect associated with the intrinsic mechanism can be 
strongly affected by carrier-impurity scattering (disorder) .L&SiiiLiLiSrfiiiiiiiSiMiiiLILiSiSiiSiiSSi (To avoid confusion, we 
use the term "intrinsic SHE" to refer to the total spin-Hall effect arising from the SO coupling terms of the Hamiltonian 
f — ■ that do not explicitly involve scattering; ultimately, this is corrected by scattering, but the part that is unaffected 
{NJ \ by scattering will be termed the intrinsic "disorder-independent" SHE.) In diffusive 2D semiconductors, there always 
exists a contribution to the intrinsic spin-Hall current which arises from spin-conserving electron-impurity scattering, 
but it is independent of impurity density within the diffusive regime. For 2D electron systems with Rashba SO coupling, 
this disorder-related spin-Hall current leads to the vanishing of the total intrinsic spin-Hall current, irrespective of the 
specific form of the scattering potential, of the collisional broadening, and of temperature^ In 2D Rashba heavy-hole 
systems, disorder affects the intrinsic SHE in a different fashion: contributions from short-range collisions to the SHE 
vanish, 18 while long-range electron- impurity scattering produces a nonvanishing disorder-related spin-Hall current, 
whose sign changes with variation of the hole densityii^tSS 
■ To date, the effect of disorder on the intrinsic spin-Hall current in p-type bulk semiconductors has been studied 
relatively little. Employing a Kubo formula, Murakami found a null disorder effect on the intrinsic SHE for short-range 
hole-impurity collisions^ The crossover of the SHE from the diffusive to the hopping regime has been investigated 
by modeling finite-size samples (with a maximum of 50 x 50 x 50 lattice sites) by Chen, et alM* In this paper, we 
^ , employ a nonequilibrium Green's function approach to study the effect of more realistic long-range hole-impurity 
scattering on the intrinsic spin-Hall current in a diffusive p-type bulk semiconductor. We find that, in such a system, 
the contribution of hole-impurity collisions to the intrinsic spin-Hall current is finite and it is independent of impurity 
density within the diffusive regime. Moreover, this disorder contribution has its sign opposite to that of the disorder- 
independent one, leading to a sign change of the total spin-Hall current as the hole density varies. Furthermore, 
we make clear that the disorder-independent spin-Hall effect arises from an interband polarization process directly 
induced by the dc electric field and it involves all hole states below the Fermi surface. In contrast to this, the disorder 
contribution to the intrinsic SHE originates from a disorder-mediated polarization between two hole bands and is 
associated only with hole states in the vicinity of the Fermi surface. Also, we numerically examine the hole-density 
dependencies of the spin-Hall conductivity and mobility. 

This paper is organized as follows. In Sec. II, we derive the kinetic equation for the nonequilibrium distribution 
function and discuss the origins of the disorder-independent and disorder-related spin-Hall currents. In Sec. Ill, 
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we perform a numerical calculation to investigate the effect of long-range hole-impurity scattering on the spin-Hall 
current. Finally, we review our results in Sec. IV. 



II. FORMALISM 
A. Kinetic equation 

It is well known that for semiconductors with diamond structure (e.g. Si, Ge) or zinc blende structure (e.g. GaAs), 
the tops of the valence bands usually are split into fourfold degenerate S — 3/2 and twofold degenerate S = 1/2 states 
due to the spin-orbit interaction (S denotes the total angular momentum of the atomic orbital). Near the top of the 
5 = 3/2 valence bands, the electronic structure can be described by a simplified Luttinger Hamiltonian 23 
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where, p = (px>Py>Pz) = (p sin P cos (f> p ,p sin 9 p sin (f> p ,p cos 9 p ) is the three-dimensional (3D) hole momentum, m is 
the free electron mass, S = (S x , S y , S z ) are the spin-3/2 matrices, and 71 and 72 are the material constants. (As in 
previous studies j&SL 2 ^^ we simplified by setting 73 = 72 in the original Luttinger Hamiltonian presented in Ref. 23). 
By a local unitary spinor transformation, U p — exp(— iS z <p p ) exp(— iS y 8 p ), Hamiltonian can be diagonalized 

as h {p) = U+h (p)U p = dia,g[e H (p),e L (p),e L (p),e H (p)]. Here, e H {p) = 71 2m 2 P 2 and £ l(p) = ^t^ P 2 are ' 
respectively, the dispersion relations of the heavy- and light-hole bands. Physically, this transformation corresponds 
to a change from a spin basis to a helicity basis. 

In a realistic 3D system, holes experience scattering by impurities. We assume that this interaction between holes 
and impurities can be characterized by an isotropic potential, V'dp — k|), which corresponds to scattering a hole from 
state p to state k. In the helicity basis, the scattering potential takes the transformed form, T(p, k) = JJ+ydp— k|)[/ k . 

We are interested in the spin-Hall current in a bulk hole system driven by a dc electric field E along the z axis. In 
Coulomb gauge, this electric field can be described by the scalar potential, V = — eE • r, with r as the hole coordinate. 
Without loss of generality, we specifically study a spin current, Jy, that is polarized along the x axis and flows along 
the y axis. In the spin basis, the conserved single-particle spin-Hall operator, jy, is defined as~ 

m-l{^p P L s x p P L +p P H s x p P H ), (2) 

with Pp and P p , respectively, as projection operators onto the states of light- and heavy-hole bands: P p = | — 
2^-(p • S) 2 , Pp = 1 — P^. Taking a statistical ensemble average, the observed net spin-Hall current is given by 

J* = £Tr[^(p)p(p)], (3) 
p 

where p(p) is the distribution function related to the nonequilibrium "lesser" Green's function, 6 < (p,cl'), as given by 
p(p) = —i J |^G < (p,o;). Also, J* can be determined in helicity basis via 

^ = £TrL?»P(p)L (4) 
p 

with jy(p) — U p jy(p)U p and p(p) = U + (p)p(p)U(p) being the helicity-basis single-particle spin current operator 
and distribution function, respectively. Explicitly, Eq. Q can be rewritten as (y6 M „ (p) are the matrix elements of p(p) 
in helicity basis; /j,, v = 1, 2, 3, 4) 

Jy = — £p{4cos 2 < ? !)pSin6lpIm[ / 5i2(p) + pu(p)] 

p 

- sin(20 p ) shi(26>p)Re[/5i 2 (p) + mip)} 
+2cos(2(/> p )cos6>pIm[ / 5i3(p) - p 2 4(p)] 

-sin(20 p )[l +cos 2 p ]Re[p 13 (p) - m(p)]} ■ (5) 
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Here, the Hermitian property of the distribution function, i.e. p(p) = p + (p), has been used. It is clear from Eq. |5| 
that contributions to the spin-Hall current arise only from those elements of the distribution function which describe 
the interband polarization, such as puip), Pi3(p): P34(p) and f>24(p)- The vanishing of spin-Hall current contributions 
from the diagonal elements of the distribution function is associated with the helicity degeneracy of the hole bands 
in p-type bulk semiconductors. The diagonal elements of the distribution function for holes in same band but with 
opposite helicities are the same, i.e. P22(p) = P33(p) and /Sn(p) = pa(p). However, the corresponding diagonal 
elements of the single-particle spin current have opposite signs due to opposite helicities, (jp22(p) = —(jy)33(p) 

and Oy)n(p) = — (jy)44(p)- As a result, the net contributions to spin-Hall current from the diagonal elements of 
distribution function are eliminated. 

In order to carry out the calculation of spin-Hall current, it is necessary to determine the hole distribution function^ 
Under homogeneous and steady-state conditions, the spin-basis distribution, p(p), obeys a kinetic equation taken the 
form, 



eE • [Vpp(p)] + *[ho(p), p(p)} = -I, (6) 



with / as a collision term given by 



i= I ^(s;,<V g p -g p - -g i:) -;s. , 7) 

G^ a,< and S p ,a,< are, respectively, the nonequilibrium Green's functions and self-energies. For brevity, hereafter, 
the argument (p,w) of these functions will be denoted by a subscript p. In the kinetic equation © above, the hole- 
impurity scattering is embedded in the self-energies, S p ' a,< . In present paper, we consider hole-impurity collisions only 
in the self-consistent Born approximation. It is widely accepted that this is sufficiently accurate to analyze transport 
properties in the diffusive regime. Accordingly, the self-energies take the forms: E p " '< = m ^ k \V{p — k)| 2 G£ a,< , 
with impurity density rii. 

It is most convenient to study the hole distribution function in the helicity basis, p(p) = J7 + (p)/5(p)J7(p), because, 
there, the unperturbed equilibrium distribution and the equilibrium lesser, retarded, and advanced Green's functions 
are all diagonal. To derive the kinetic equation for the helicity-basis distribution, p(p), we multiply Eq. © from 
left by U£ and from right by U p . Due to the unitarity of U p , the collision term in the helicity basis, /, has a form 
similar to Eq. J7J, but with the helicity-basis Green's functions and self-energies, G p ' a,< = J7 + (p)G p a,< £/(p) and 
£r,a,< _ [/+(p)xi p ^<(7(p) 5 respectively, replacing those of the spin-basis, G p ,a,< and S p ,a,< . The left hand side 
(LHS) of Eq. © is simplified by using the following facts: f/+V p/ 5(p)C/ p = V p/ 5(p) - V p C7+f7 p/ 5(p) - p(p)U+V p U p 
and VpUpUp — —U p W p U p . Thus, the kinetic equation in helicity basis may be written as 

eE - {V p/ 5(p) + [p(p),V P U+U p ]}+i[h (p),p(p)] = -I. (8) 

In this equation, the helicity-basis self-energies, S p ,a,< , take the forms, 

Sp' a,< =^^f(p,k)G^ < f+(p,k). (9) 

k 

In this paper, we restrict our considerations to the linear response regime. In connection with this, all the functions, 
such as the nonequilibrium Green's functions, self-energies and distribution, can be expressed as sums of two terms: 
A = Aq + Ai, with A as the Green's functions, self-energies or distribution function. Aq and A\, respectively, are the 
unperturbed part and the linear electric field part of A. In this way, the kinetic equation for the linear electric field 
part of the distribution, pi (p) , can be written as 

eE • V p/ 5 (p) - eE • [p (p), U+V P U P ] + i[h (p), pi(p)] = -/« (10) 

with as the linear electric field part of the collision term /: 

/M = 



2^ 



^lp^Op "+" ^lp^Op — '-'lp^Op — '-'lp^Op 



+^0pGi p + Sg p G° p — GgpSfp — G^ p S"p . (11) 

Further, we employ a two-band generalized Kadanoff-Baym ansatz (GKBA)2£i2I to simplify Eq. I|1U|) . This ansatz, 
which expresses the lesser Green's function through the Wigner distribution function, has been proven sufficiently 



4 



accurate to analyze transport and optical properties in semiconductors^ To first order in the dc field strength, the 
GKBA reads, 

Gfp = -GSpPi(p) + Pi(p)G a p - G^ p po(p) + A)(P)G? P , (12) 

where the equilibrium distribution, and retarded and advanced Green's functions are all diagonal matrices: po(p) = 
d\ag[n F (eH{p)),n F (e L (p)),n F (£ L (p)),n F (e H (p))} and Gg' a (p) = diag[(w - e H {p) ± iS)' 1 , (w - e L (p) ± i<5) _1 ,(w - 
Sl(p) ± «<5)~\ (w — £h(p) ± w ith the Fermi function np(w). We note that G^ p in the collision term leads to a 

collisional broadening of the nonequilibrium distribution. In the present transport study, such collisional broadening 
plays a secondary role and can be ignored. Based on this, the collision term, JW, no longer involves the linear electric 
field part of the retarded and advanced Green's functions. 

It is obvious that the driving force in Eq. (|lt)(l comprises two components: the first of which, eE • V p po, is diagonal, 
while another one, — eE • [po(p), E^VpZ7 p ], has null diagonal elements. In connection with this, we formally split the 
kinetic equation into two equations with pi(p) — p{(p) + Pi 1 (p) as 

eE ■ V p/ 3 (p) +i[h (p),p{(p)] = (13) 



- eE • [p (p), U+ V p C/ p ] + i[fto(p), pVip)} = 0, (14) 

wherein p{(p) and /5{ 7 (p) can be approximately determined independently, as discussed below. We note that the 
solution of Eq. 114(1. p 1 ^ (p), is off-diagonal and independent of impurity scattering. The matrix elements of p{ (p) will 
be denoted by (p{' 77 ) Ai;y (p), and from Eqs. and (JSJ, we correspondingly write spin-Hall conductivity contributions 
based on J* = J* | J + J* \ JI as 

(°% = J?! 7 IE = J2 Tr[^(p)/5((p)]; 
p 

<p JI )i,= /E = j2nfy(p)pV(p)}- (15) 
p 

It is evident that the diagonal driving term of Eq. I|13|) . eE • V p /5o, is free of impurity scattering. Since [ho, p{(p)] 
is off-diagonal, the diagonal parts of this equation lead to diagonal p[{p) elements, (/S{) MM (p) (p — 1...4), of order of 
(rii) -1 in the impurity density. Substituting these diagonal elements, (p{) ^(p) , into the off-diagonal elements of the 
scattering term, I^, and considering the fact that the terms on LHS of the off-diagonal components of Eq. I|13|) are 
proportional to the off-diagonal elements of p[ (p), we find that the leading order of the off-diagonal elements of p[(p) 
in the impurity-density expansion is of order (rii) , i.e. independent of rij. This result implies that, in general, there 
always exists a contribution to the spin-Hall current which is disorder-related but independent of impurity density 
within the diffusive regime. On the other hand, the off-diagonal impurity-density-independent p[(p) elements, as 
well as all the nonvanishing elements of /5{ 7 (p), make contributions to the scattering term, /W, which are linear in 
the impurity density, while the terms involving diagonal elements, (p{) w (p), are independent of nj. Hence, the 
contributions to 1^ from off-diagonal elements of pi(p) can be ignored and effectively involves only the diagonal 
elements of the distribution. Correspondingly, Eqs. I|13J) and (|14|l are approximately independent of each other and 
can be solved separately. 



B. Disorder-independent spin-Hall effect 

The disorder- independent spin-Hall current is associated with p{ 7 (p), the solution of Eq. I(14J1 . The nonvanishing 
elements of this function are given by 

G3f)l2(p) = -(p( / )2l(p) = (p( / )34(p) = -(/5( / ) 43 (p) 

= ^ le Esm9 p [f H (p)-f a L (p)}, (16) 

with /(^(p) = n F [eH(p)] and (p) = n F [eL{p)\, while its remaining elements, such as (/5{ 7 )i3(p), (/5{ / )24(p), etc. 
vanish. Substituting p{ ! (p) into Eq. JSJ, we find that the disorder- independent contribution to intrinsic spin-Hall 
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current, J? , can be written as 

771 />00 

J y\ U = -^2j o [fo(p)-fo(p)}dP' (17) 

This result agrees with that obtained in Ref.Eil 

Obviously, the nonvanishing of J*\ is associated with the nonzero driving term on LHS of Eq. 114|l . which is just 
the interband electric dipole moment between the heavy- and light-hole bands. Thus, the disorder-independent spin- 
Hall effect arises essentially from the polarization process between two hole bands directly induced by the dc electric 
field. Such a polarization can also be interpreted as a two-band quantum interference process. It should be noted that 
this polarization process affects only those off-diagonal p{ ! (p) elements which describe dc- field induced transitions 
between hole states in the light- and heavy-hole bands. Of course, such transition processes are not restricted only to 
hole states in the vicinity of the Fermi surface: they contribute from all the hole states below the Fermi surface. As a 
result, the disorder-independent spin-Hall current given by Eq. ljT7|l is a function of the entire unperturbed equilibrium 
distribution, np(w), not just of its derivative, dn-p(u)/duj, at the Fermi surface. 



C. disorder- related spin-Hall effect 

To simplify Eq. (|I3|I . we first analyze symmetry relations between the elements of the distribution function p{(p) in 
the self-consistent Born approximation. Since the distribution function is a Hermitian matrix, only the independent 
elements (/5{) Atl /(p) with p, v = 1...4 and p <v need to be considered. We know that (/5{)ii(p) and {p\)^{p) describe 
the distributions of the heavy holes having spins S z — 3/2 and S z = —3/2, respectively. In equilibrium, heavy hole 
populations in degenerate states with S z = 3/2 and S z = —3/2 distribute equally. Out of equilibrium, the dc electric 
field action on these hole populations is also the same. Hence, the noncquilibrium distribution of the heavy holes 
with S z — 3/2 is the same as that of the heavy holes with S z — —3/2, i.e. (p[)n(p) — (^{^(p). An analogous 
relation for light holes is also expected to be valid: (p I 1 )22(p) — {p{)33(p)- Indeed, substituting these symmetrically 
related diagonal elements of the distribution p[(p) into the scattering term, we find l[i — 1^ , — and 

^23 } = I32 = = = °> which 

are consistent with the elements on the LHS of Eq. 11311 . As another consequence 
of these relations ((/5{)n(p) = (/5{) 44 (p) and (/5{) 22 (p) = (Pi)33(p)), we also obtain symmetry relations between the 
remaining off-diagonal elements of J^: = — ^34 and = 7^4 , which result in symmetry relations for the p[(p) 
elements as: (/5{)i 2 (p) = {p{)34(p) and (p[)i 3 (p) = — ( ( o{) 2 4(p). Hence, to determine the disorder-related spin-Hall 
effect, one only needs to evaluate the diagonal elements, (p[)n(p) and (p{)22(p), and the off-diagonal elements, 
(Pi)i2(p) and (pf)i 3 (p). 

From Eq. (|13|) , it follows that the diagonal p\ (p) elements are determined by the integral equation 

-eE- V p n F [e^p)} = vr^ \V(p - k)| 2 { ai (p, k)[(p() w (p) 

k 

-(p[) w (k)]A w + a 2 (p,k)[(p[) w (p) 

-(p[) w (k)A^]. (18) 

Here, fx = 1, 2, respectively, correspond to the heavy- and light-hole bands: £i{p) = £h(p), £2(p) = £l(p), ft = 3 — fi, 
= 5[e^(p) — e v (k)]. The factors ai(p, k) and a2(p, k) are associated only with the momentum angles: 

ai(p,k) = i{2 + 6 cos 2 </> p fc[sin 2 9 p - cos 2 k ] 
+6 cos 2 9 P cos 2 6*k[l + cos 2 4> p k] 

+'3cos(f) pk cos(29 p )cos(29 k )}, (19) 



a 2 (p,k) =2-a 1 (p,k), (20) 

where tfi p k = <fi P — 0k- From Eq. i|18|) . we see that we may remove the dependence of (Pi)/^(p) on momentum angle 
4>p by redefining the angular integration variable as 0k — > fipk = 0p — 0k, taken jointly with the facts that the left 
hand side does not depend on </> p and the potential V{p — k), as well as the factors ai(p, k) and a 2 (p, k), depends on 
p and 0k only through the combination . 

Analyzing the components of the scattering term in the kinetic equation for the off-diagonal elements, (p{)i 2 (p) 
and (/5{)i3(p), we find that these elements of the distribution p{(p) are similarly effectively independent of p . In 
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connection with this, contributions to the disorder-related spin-Hall current, J* \ T , from (pf)i3(p) an d R e [(Pi)i2(p)] 
vanish under the </> p -integration in Eq. J^J), and only the imaginary part of (p{)i2(p) makes a nonvanishing contribution 
to Jy\ ■ Hence, 



JZ\* = ^^Ep{ cos2 ^ sin0 P Im [^)i2(P)]} . ( 21 ) 
p 

with 

lm[(p[) 12 (p)] = ^ E in P -k)| 2 a 3 (p,k) 

~ 2P k,ft=l,2 

x(-l)"{A w [(pi) w (p)-(pf) w (k)] 

-A^[(p() w (p) - (pi) w (k)]}, (22) 

and 

a 3 (p,k) = -i{sin(26' p )[cos 2 6' k -sin 2 6' k cos 2 (/)p fe ] 

+ sin(26» k )cos0 pfc [l - 2 cos 2 p ]}. (23) 

From Eqs. (|18|l and 1)22(1 . we see that Jy ^ is independent of impurity density. In contrast to the disorder- independent 
case, the disorder-related spin-Hall current involves only the derivative of the equilibrium distribution function, i.e. 

dnF(u>)/duj. This implies that Jy | is constituted of contributions arising only from hole states in the vicinity of the 
Fermi surface, or in other words, from hole states involved in longitudinal transport. Physically, the holes participating 
in transport experience impurity scattering, producing diagonal p{(p) elements of order of n~ x . Moreover, the 
scattering of these perturbed holes by impurities also gives rise to an interband polarization, which no longer depends 
on impurity density within the diffusive regime. It is obvious that in such a polarization process the disorder plays 
only an intermediate role. It should be noted that Jy | generally depends on the form of the hole-impurity scattering 
potential, notwithstanding its independence of impurity density in the diffusive regime. 

The fact that the total spin- Hall current, Jy = J* | + J®\ , consists of two parts associated with hole states 
below and near the Fermi surface, respectively, is similar to the well-known result of Stfeda 2 - in the context of the 2D 
charge Hall effect. In 2D electron systems in a normal magnetic field, the off-diagonal conductivity usually arises from 
two terms, one of which is due to electron states near the Fermi energy and the other is related to the contribution 
of all occupied electron states below the Fermi energy. A similar picture has also recently emerged in studies of the 
anomalous Hall effect ,22i2i 



III. RESULTS AND DISCUSSIONS 

To compare our results with the short-range result presented in Ref. l2lL we first consider a short-range hole- impurity 
scattering potential described by: V(p— k) = u, with u as a constant. Substituting Eq. i|22|) into Eq. J5J) and performing 
integrations with respect to the angles of p or k, respectively, for terms involving (/5{) MA1 (k) or (/5{) MM (p), we find 

that the contribution of short-range disorder to the spin-Hall current vanishes, i.e. J*\* = 0. This implies that for 

short-range hole- impurity collisions, the total spin-Hall current is just the disorder- independent one, Jy — Jy^ 1 ■ 
This result agrees with that obtained in Ref.l2ll 

Furthermore, we perform a numerical calculation to investigate the effect of long-range hole-impurity collisions 
on the spin-Hall current in a GaAs bulk semiconductor. The long-range scattering is described by a screened 
Coulombic impurity potential V(p): V(p) — e 2 / (eoe)[p 2 + with e as a static dielectric constant^ do is 

a Thomas-Fermi-Debye type screening length: d 2 D = it 2 e e / (e 2 y / 2m 3 E F )2^ 1 /' i [{^i + 2j 2 )~ 3 ^ 2 + (li — 272)~ 3 / 2 ]~ 2,/3 , 
with Ep ~ (37T 2 A r p /2) 2 / 3 /(2m). The material parameters 71 and 72 are chosen to be 6.85 and 2.5, respectively^ In 
our calculation, the momentum integration is computed by the Gauss-Legendre scheme. 

In the present paper, we address the spin-Hall effect at zero temperature, T = 0. In this case, the disorder- 
independent spin-Hall current can be obtained analytically from Eq. ((17(1 : J*\ = eE[kp — kp]/(6Tr 2 ), with kp and 
kp as the Fermi momenta for heavy- and light- hole bands, respectively. In order to investigate the disorder-related 
spin-Hall effect, we need to compute the distribution function p{(p) at the Fermi surface. In this calculation, we 
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FIG. 1: Hole-density dependencies of (a) total Oy Z and fiy Z , and (b) disorder-independent (a !I )y Z and {p II )y Z , in a bulk GaAs 
semiconductor. The material parameters for GaAs are: 71 = 6.85 and 72 = 2.5. The lattice temperature is T = OK. 

employ a "singular value decomposition" method 34 to solve the integral equation, Ea. (|18|l . for the diagonal p[{p) 
elements. The obtained diagonal elements are then employed to determine Im[( j o{)i2(p)] using Eq. i|22|l . Following 
that, we obtain the disorder-related spin-Hall current from Eq. (|21|l . performing the momentum integration. 

In Fig. 1, the calculated total and disorder-independent spin-Hall conductivities, a% z = Jy/E and ((J n )y Z = 

Jy\ I! /E, and the total and disorder-independent spin-Hall mobilities, [1* — (jy Z /N p and (fJ- 11 )^ = (a n )y Z /N p , 
are shown as functions of the hole density. The spin-Hall mobility, analogous to the mobility of charge transport, 
characterizes the average mobile ability of a single spin driven by the external field. This quantity has the same units 
in 2D and 3D systems. 

From Fig. 1, we see that, with increasing hole density, the total spin-Hall conductivity first increases and then 
decreases and even becomes negative as the hole density becomes larger than N pc = 3 x 10 24 m~ 3 . This behavior 
of the hole-density dependence of total spin-Hall conductivity is the result of a competition between the disorder- 
independent and disorder-related processes. The contributions to spin-Hall conductivity from these two processes 
always have opposite signs and their absolute values increase with increasing hole density. Considering total spin-Hall 
conductivity, the disorder-related part, (c 7 )y Z , is dominant for high hole density, while (<r 77 )^ z is important in the 
low hole-density regime. Notwithstanding this hole-density dependence of c^ 2 , the total spin-Hall mobility, fj,y Z , as 
well as the disorder-independent one, monotonically decreases with increasing hole density. 

It should be noted that the total spin-Hall mobility in bulk systems has the same order of magnitude as that in 2D 
hole systems. We know that the spin-Hall conductivity in 2D hole systems is of order of e/iri& For a typical 2D hole 

density, n p 2D ^ = 1 x I0 12 cm~ 2 , the corresponding spin-Hall mobility is about 0.05m 2 /Vs. 

In the present paper, we have ignored the effect of collisional broadening on spin-Hall current. Since Jy\* is 
associated only with the hole states in the vicinity of the Fermi surface, the neglect of broadening in the disorder- 
related spin-Hall current is valid for eft > 1 (sf is the Fermi energy and r is the larger of the relaxation times 
for holes in the different bands at the Fermi surface, tl(sf) and th(sf)- t — max[rx(£i?), th (sf)])- This condition 
coincides with the usual restriction on transport in the diffusive regime and is satisfied for p-type bulk GaAs with 
mobility approximately larger than lm 2 /Vs (for N p > 5 x 10 22 m~ 3 ). On the other hand, the disorder-independent 
spin-Hall conductivity involves contributions from all hole states in the Fermi sea and hence it may be strongly 
affected by collisional broadening. To estimate the broadening effect on the disorder-independent SHE, we add an 
imaginary part to ho(p) and use ho(p) +«7(p) instead of ho(p) in Eq. lfH|) ( j(p) is a diagonal matrix describing the 
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broadening: (7)11 (p) = {^uip) = l/2T H (e H {p)) and (7) 22 (p) = (7)33(p) = 1/ 2t l( £ £ (??)))• In this way, J*\ takes 
a form similar to Eq. {T7|) but with an additional factor, (2'y 2 p 2 ) 2 / {[^P 2 ] 2 + [1 / '2t h (e h (p)) - 1/2tl(£l(p))] 2 }, in the 
momentum integrand. Performing a numerical calculation, we find that, in the studied regime of hole density, the 
effect of collisional broadening on the disorder-independent spin-Hall current is less than 1% for p-type bulk GaAs 
semiconductors with mobility approximately larger than 5 m 2 /Vs. Thus, in such systems, the effect of collisional 
broadening on the total spin-Hall conductivity can be ignored. It should be noted that in our calculations, we 
computed tl } h (^) by considering short-range hole-impurity scattering: 1/tl,h(s) = 2imiU 2 VL,H{s) with the densities 
of hole states in the light- and heavy-hole bands taken as vl,h( £ ) = 2^2 p 5(s — sl,h(p))- The quantity rnu 2 is 
determined from the mobility of the system: [i = e[N^T[ J {£F)/mi J + th{£f) /™h\/N p , where — m/^i + 272) 
and m H = m/(jx - 272) are the effective masses of holes and N£ /iVjf = [(71 - 27 2 )/(7i + 27 2 )] 3 / 2 with and 
being the hole densities in the light- and heavy-hole bands, respectively. 

On the other hand, in our considerations, the impurities are taken to be so dense that we can use a statistical 
average over the impurity configuration. This requires that < L (L is the characteristic size of the sample and 
Lb is the larger of the diffusion lengths of holes in the light- and heavy-hole bands). Failing this, the behavior of the 
holes would become ballistic, with transport properties depending on the specific impurity configuration. 

IV. CONCLUSIONS 

We have employed a nonequilibrium Green's function kinetic equation approach to investigate disorder effects on the 
spin-Hall current in the diffusive regime in p-type bulk Luttinger semiconductors. Long-range hole- impurity scattering 
has been considered within the framework of the self-consistent Born approximation. We have found that, in contrast 
to the null effect of short-range disorder on the spin-Hall current, long-range scattering produces a nonvanishing 
contribution to the spin-Hall current, independent of impurity density in the diffusive regime. This contribution has 
its sign opposite to that of the disorder-independent one, leading to a sign change of the total spin-Hall current as 
the hole density varies. We also made clear that the disorder-independent spin-Hall effect arises from a dc-field- 
induced polarization associated with all hole states in the Fermi sea, while the disorder-related one is produced by a 
disorder-mediated polarization and relates to only those hole states in the vicinity of the Fermi surface. The numerical 
calculation indicates that with increasing hole density, the total spin-Hall mobility monotonically decreases, whereas 
the spin-Hall conductivity first increases and then falls. 

In addition to Jy, we also examined other components of the spin current. We found that the previously discovered 
"basic spintronics relation"^ which relates the ith component of the spin current along the direction j, Jj, and the 
applied electric field, Ek, by Jj = <r s £ijkEk with etjk as a totally antisymmetric tensor, still holds in the presence of 
spin-conserving hole-impurity scattering. 
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